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Abstract. We suggest and study the stable disk- and cigar-shaped gap solitons of 
a dipolar Bose-Einstein condensate of ^^Cr atoms localized in the lowest band gap by 
three optical-lattice (OL) potentials along orthogonal directions. The one-dimensional 
version of these solitons of experimental interest confined by an OL along the dipole 
moment direction and harmonic traps in transverse directions is also considered. 
Important dynamics of (i) breathing oscillation of a gap soliton upon perturbation and 
(ii) dragging of a gap soliton by a moving lattice along axial z direction demonstrates 
the stability of gap solitons. A movie clip of dragging of three-dimensional gap soliton 
is included. 
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Solitons in Bose-Einstein condensate (BEC) have drawn a great deal of attention 
recently. Experimentally, bright solitons were created in attractive BEC of ^Li [T] and 
^^Rb [2] atoms resulting from a cancellation of nonlinear attraction and linear dispersive 
effects. For a repulsive BEC, dark solitons have also been experimentally observed [3]. In 
a repulsive BEC on a periodic optical-lattice (OL) locahzed (dynamically) stable states 
of gap solitons in the band gap can be made [HE]. An OL simulates the periodic electron- 
atom potential in a solid, and the study of gap sohtons in OL and the related band-gap 
structure in the presence of nonlinear interaction is also of interest in condensed-matter 
physics [6]. The mathematical model of gap solitons and atomic band-gap in a BEC 
is identical to that of light propagation in photonic crystals with cubic nonlinearity in 
presence of photonic band-gap and this makes the study of gap solitons in a BEC of 
general interest [1] . Later, the possibility of gap solitons in Fermi superfluid has been 
suggested [7]. Gap solitons have been observed in photonic crystals [8]. A gap soliton 
of about 250 ^^Rb atoms in a nearly one-dimensional (ID) OL (Qj [10] has been created 
experimentally. In the experiment, the OL was subjected to acceleration, in order to 
push atoms into a gap soliton state [9] . Other possibihty for the creation of gap solitons 
may employ the phase-imprinting method 

The alkali metal atoms used in BEC experiments have negligible dipole moment. 
However, most bosonic atoms and molecules have large dipole moments and a ^^Cr BEC 
[12^ . with a larger long-range dipolar interaction superposed on the short-range atomic 
interaction, has been realized. Thus one can study a dipolar BEC (DBEC) with variable 
short-range interaction using a Feshbach resonance [I3]. The stability of a DBEC 
depends not only on the atomic interaction, but also on trap geometry [El [HI [15]. In a 
disk configuration the dipolar interaction is repulsive and a DBEC is more stable; while, 
in a cigar configuration the dipolar interaction is attractive and a DBEC is less stable 
due to collapse instability [T^ [TU] . The controllable short-range interaction together 
with the dipolar interaction makes the DBEC an attractive system for experimental 
soliton generation and a challenging system for theoretical investigation. There have 
already been studies of bright solitons in DBEC [T7] . 

Using the numerical and variational solution of the CP equation we predict 
and study stationary and dynamical properties of three-dimensional (3D) DBEC gap 
solitons, of both cigar and disk shapes, localized in the lowest band gap by three 
orthogonal OL potentials. The BEC gap solitons are realizable only for all repulsive 
atomic interactions below a maximum value so that the chemical potential falls in the 
band gap. The DBEC gap solitons in disk shape can also be formed for weakly attractive 
short-range interaction. The cigar-shaped DBEC gap solitons are formed only for short- 
range repulsion above a limiting value. We also consider the one-dimensional gap solitons 
of direct experimental interest confined by an OL along the dipole moment direction 
and harmonic traps in transverse directions. 

We numerically explore the breathing oscillation of the gap solitons upon a small 
perturbation. We illustrate the dragging of the gap solitons by an OL moving along 
the axial z direction [18], [19]. These solitons can be dragged without deformation for a 
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reasonably large velocity. 

Here we study the DBEC gap soliton of atoms, each of mass m, using the 
dimensionless GP equation |12] 

+ N j [/rf,(r-r')|0(r',t)|W 

with dipolar interaction Uddi^) = 30^^(1 — 3cos^6') / , R = r — r'. Here V^(r) 
is the confining potential, 0(r, t) the wave function at time t with normalization 
/|0(r,t)p(ir = 1, a the atomic scattering length, 6 the angle between R and the 
polarization direction z. The constant add = ^qJj? m / {1271^^) is a length characterizing 
the strength of dipolar interaction and its experimental value for ^^Cr is 15ao [12], 
with oq the Bohr radius, Ji the (magnetic) dipole moment of a single atom, and /ig 
the permeability of free space. The OL potential in a specific direction, say is 
V{z) = s zE ji sin^ {2tt z / X), where E^i = /i^/(2mA^) is the recoil energy of an atom, 
A is the wave length of the laser and is the strength of the OL. In ([T]) length is 
measured in units of A/27r (taken here as 1 fim for a far infrared laser), time t in units 
of mA^/27r/i, and energy in units of /i^/mA^. The dimensionless 3D periodic OL trap 
can now be written as 

V^(r) = -[Vp{cos (2x) + cos {2y)} + cos {2z)], (2) 

where the parameters Vp and Vz are the strengths of the OL's in radial and axial 
directions and can be varied to achieve the disk- {Vz > Vp) and cigar-shaped (Vp > Vz) 
gap solitons in DBEC. 

The gap solitons considered here predominantly have a Gaussian shape and for 
these solitons a Gaussian variational solution is known to lead to a good description 
[7] . The variational approximation provides an analytical understanding and also yields 
interesting results when the numerical procedure is difficult to implement. A better 
but more complicated description can be obtained with a different ansatz of variational 
functions [20]. The Lagrangian density of ([T]) is given by 

C= {<P<Pt - + ^iV^r + V{v)\<P\' + 27raAr|0|4 
N r 

+ Udd{r-r')\<p{r')\'dr'. (3) 

We use the Gaussian ansatz [2T] 

7r-3/4 / 2 ^2 \ 

0(r, t) = — exp - — + lap^ + t^z^ (4) 

for a variational calculation, where Wp and Wz are time-dependent radial and axial 
widths, and a and (3 time-dependent phases. The effective Lagrangian L = J Cdr (per 
particle) becomes 

L = {w^a + wlP/2) - 2Vpexp{-wl) - Vz exp{-wl) 
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+ l/{2wl) + l/{Awl) + 2wy + wl(3^ 

+ Nadd/{V2mwlw^)[a/add- f{K)], (5) 



with 



/(k) = [l + 2/t2-3/t2rf(/t)]/(l-K^), (6) 

din) = == — , K= (7) 

The Euler-Lagrange equations for parameters Wp.Wz.cy. and [3 yield the following 
equations for widths Wp and Wz 

^^^^[2a-a^ (8) 
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with 

(^(ft:) = [2 - 7k' - + 9K^d{K)]/{l - K^f, (10) 
h{K) = [1 + lOfi:' - 2k^ - 9K^d{K)]/{l - K'f. (11) 

The chemical potential /i is given by 

^1^1^ 2N[a-addf{n)] 
^ 2wf, Awl V2TrWzWj 

- 2Vpexp{-wl) - Vzexp{-wl). (12) 

For gap solitons the system must be repulsive. For a normal BEC {add = 0), 
attraction corresponds to a < and repulsion to a > and gap solitons are formed 
below a limiting repulsive scattering length < a < ac- For add > 0, in the disk shape, 
the dipole moment contributes repulsively and due to this extra repulsion a DBEC gap 
soliton can be formed in a window of scattering lengths —ai < a < a2 between a limiting 
attractive (— ai) and repulsive (02) limits. However, in a weak cigar shape, the dipole 
moment contributes attractively and due to the extra attraction a DBEC gap soliton 
can be formed in a window of repulsive scattering lengths 03 < a < 04. The limiting 
values ai, 02, as and a^ can be obtained from a solution of the GP equation. In a strong 
cigar regime, the dipole moment contributes to strong attraction and no gap soliton can 
be formed due to collapse instability. 

The cigar-shaped quasi-lD gap solitons with a radial harmonic trap of frequency 
fip satisfy ([T ) with V{r) = f2^p'/2— cos(22;). In this case it is convenient to solve the 
ID GP equation [221 123] with the reduced dipolar interaction 

U}^iZ) = 3addN[A6iVi)/3 + 2^1- v^(l + 2t) 

xe*{l-Erf(y^)}]/(v^rfp)^ (13) 

where Z = \z — z'\, dp = 1/^JVt'p, t = \Z j [y/2dp)f and where we have included the 
proper 5-function term mentioned in p2]. The ID GP equation is given by 
.d4>xD{z,t) r 1 92 2aAr 2 
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Figure 1. (Color online) Band (shaded area), gap (white area) and the numerical (N) 
and variational (V) a — fi plots of 3D DBEC gap solitons for add = 15ao, = 500 and 
Vp = 5, (a) 14 = 50 (disk) and (b) = A (cigar). 
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Using the Gaussian ansatz 
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the variational Lagrangian becomes 
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The variational equation for Wz becomes 

,3 



idz + 



The corresponding chemical potential is 



1 _^ 2N[a - addh{Ko)] 
wi \/2iw1dl 



fx 
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(14) 



(15) 



(16) 



(17) 



4^2 



27nUzd'i 



exp(M;2) ■ 

To study the gap solitons, we perform a 3D numerical simulation employing 
real-time propagation with Crank-Nicolson method [23]. The numerical solutions are 
obtained by averaging the wave function over iterations with the variational solution as 
input. The dipolar interaction is evaluated by fast Fourier transform [H] . 

The stationary gap solitons appear in the band gap of the OL and a prior knowledge 
of the band and gap of the periodic OL is of advantage. The dipolar interaction is most 
prominent in the cigar (attractive) and disk (repulsive) shapes and hence we study the 
gap solitons in these two shapes. The OL parameters in these cases are taken as (a) 
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Ip = 5, \4 = 50 (disk shape) and (b) = 5, Vz = 4 (cigar shape). The band and 
gap in these cases are shown in figure [T] together with the variational and numerical 
chemical potential /x for the stable gap solitons in the lowest band gap as a function of 
scattering length a for a DBEC {add = 15ao). The gap solitons in the higher band gaps 
are found to be dynamically unstable. In figure [T| the agreement between variational 
and numerical results worsens for large a. The system is more repulsive for large a and 
tends to occupy multiple OL sites and a single-peak Gaussian variational ansatz may 
not provide a good approximation to reality. In the cigar shape {Vz = 4), the dipolar 
interaction is attractive and hence DBEC gap solitons are possible for a = > 0.1343 
nm as in figure [l] (b). In the disk shape {Vz = 50), the dipolar interaction is repulsive, 
and hence DBEC gap solitons are possible for a = —ai > —0.2635 nm as in figure [T] 
(a). 

In figure |2] we show the 3D contour density profiles of some typical stationary gap 
solitons with = 500, Vp = 5, add = ISoq for variable a and V^. We show disk-shaped 
profiles for Vz = 50 and (a) a = 0.5 and (b) a = —0.1 and cigar-shaped profiles for Vz = 4: 
and (c) a = 1 and (d) a = 0.5. For both disk and cigar shapes, the increased scattering 
length implies more nonlinear repulsion and one can have more atoms in adjacent OL 
sites. Although the gap solitons in figure [2] occupy multiple OL sites, more than 95 % 
of the matter is contained in the central site and an approximate Gaussian distribution 




Figure 2. (Color online) 3D contour density plot of gap solitons for ~ 500, Vp = 
5, add = 15ao. (a) K = 50,a = 0.5, (disk), (b) Vz = 50, a = -0.1, (disk), (c) 
= 4,a = 1, (cigar), (d) Vz — 4:,a = 0.5, (cigar). The density at contour is 0.002. 
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Figure 3. (Color online) (a) Band (shaded area), gap (white area) and the numerical 
(N) and variational (V) a—// plot of ID DBEC gap solitons for Q,p = l (disk) and Vlp = A 
(cigar). Variational and numerical densities of the ID DBEC for (b) Op = 4, a = 0.75 
nm, (c) rip = 4, a = 2.75 nm, (d) Jlp = 1, a = 1 nm, and (e) ilp = 1, a = 3 nm. In all 
cases add — 15ao, N — 500 and Vz — 5. 



is valid. 

Next we study ID gap solitons trapped by the OL —Vz cos(22;) along the z direction 
and harmonic trap ^^pP^ in transverse directions as calculated from numerical and 



variational solutions of the reduced ID GP equation [22] with dipolar interaction ( 13 ). In 



figure |3] (a) we show the bands and gaps as well as the a — yU plot of the ID gap solitons. 
For fip = 1, the DBEC is of disk shape and the contribution of the dipolar term is 
repulsive, and for Vtp = 4, DBEC is of cigar shape and the dipolar term contributes 
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Figure 4. (Color online) Numerical (num) and variational (var) rms sizes {p),{z) 
during breathing oscillation of a disk-shaped 3D DBEC for N = 500, add = ISooj = 
50, Vp = 5 as a is changed from to 0.1. 
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Figure 5. (Color online) Snapshots of 3D contour density profiles of a cigar-shaped gap 
soliton for Vj — 4, = 5, = 500, add — 15ao, and a — 0.5 nm during dragging by the 
moving OL -V^p[cos(2a;)+cos(2?/)]-2l4 cos[2{z-vt)],v = 0.75 at times i = 0, 4, 8, 12, 16 
and 20. 

attractively. In figures [s] (b) — (e) we plot the typical numerical ID density profiles of 
some gap solitons and compare with variational results. Although, variational results 
exist across the bands, stable DBEC 3D gap solitons exist in the lowest band gap away 
from the band. In ID some stable gap solitons are found in the first excited gap. Gap 
solitons cannot be stabilized close to the bands, e. g., near /i = —45 and —49 in figure 
[l] (a), near /i = —4.5 and —8 in figure [l] (b), and near /z = 1 and 4 in figure [sj In figures 
[l](a) and [3] (a) we clearly see that gap solitons are possible for small negative scattering 
lengths in disk-shaped DBEC and that, in the second band gap, ID gap solitons can be 
stabilized only in a small domain of fi values between two bands. 

To test the stability of the 3D DBEC gap solitons we now consider the breathing 
oscillation of a disk-shaped soliton initiated by slightly changing the scattering length 
a. Such a change in a can be made by varying a background magnetic field near 
a Feshbach resonance [I3] . In figure |4] we plot the variational and numerical rms sizes 
(p), (z) versus time t. We find, using Fourier analysis, that the principal axial frequencies 
1.0137 (variational) and 1.0812 (numerical) and the principal radial frequencies 0.2421 
(variational) and 0.2469 (numerical) are in good agreement with each other. 

Finally, we study the stability of a cigar-shaped 3D DBEC when dragged by an 
OL moving in the axial z direction. It has been experimentally found that a BEC 
confined by a transverse harmonic trap remains stable |19] when dragged by a moving 
ID OL along the axial z direction below a critical velocity Vc of half the recoil velocity 
[vr = h/{m\)], which in present units (A = 27r and m = h = 1) is Vc = vr/2 = 0.5. 
Similar result is found to be true in the case of present 3D DBEC gap solitons. Steady 
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dragging is possible in this case even for velocities slightly larger than this limit. Instant 
snapshots of dragging dynamics with a velocity v = 0.75(< Vc = 1) along z direction is 
shown in figure [5] where we find that the gap soliton can be dragged without distortion 
with a reasonably large velocity for a considerable time. A movie of dragging of 
soliton of figure [5] showing stability is prepared and contained in supplementary clip SI 
(see supplementary data available at stacks.iop.org/JPhysB/44/xxxxxx/mmedia, also 
available at http://www.youtube.com/watch?v=wJbbKmsjQuU). 

To conclude, we suggested the possibility of 3D DBEC gap sohtons of about 1000 
Cr atoms confined in the lowest band gap by three OL in orthogonal directions and 
studied their statics (shape and chemical potential) and dynamics (breathing oscillation 
and dragging by an OL). In addition we studied ID DBEC gap solitons using reduced 
ID GP equations with a transverse harmonic trap and an axial OL along polarization 
direction. The 3D DBEC gap solitons are found to be dynamically stable during 
breathing oscillation and dragging for a long enough time for experiments and, with 
available technology, these solitons could be created and studied in laboratory. The 
present study opens up new directions of research that include, among others, excited 
states of gap solitons |1] and vortex gap solitons [2S] in DBEC. 
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